Ab initio study of the Ti02 Rutile(110)/Fe interface 
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Adsorption of Fe on the rutile (llO)-surface is investigated by means of ah initio density functional 
theory calculations. We discuss the deposition of single Fe atoms, an increasing Fe coverage, as well 
as the adsorption of small Fe clusters. It is shown that the different interface structures found in 
experiment can be understood in terms of the adsorption of the Fe atoms landing first on the rutile 
surface. On the one hand, strong interface bonds form if single Fe atoms are deposited. On the other 
hand, the Fe-Fe bonds in deposited Fe clusters lead to a three-dimensional growth mode. Mainly 
ionic Fe-oxide bonds are formed in both cases and the electronic band gap of the surface is reduced 
due to interface states. Besides the structural and electronic properties, we discuss the influence of 
the interface on the magnetic properties finding stable Fe moments and induced moments within 
the interface which leads to a large spin polarization of the Fe atoms at the rutile (110)/Fe interface. 



I. INTRODUCTION 

The (llO)-surface of rutile is one of the most im- 
portant metal-oxide surfaces and has been subject of a 
large number of theoretical and experimental work in the 
last decades.^ ^ Numerous applications ranging from thin 
film capacitors to optical devices are based on the excep- 
tional dielectric properties of this oxidic- surf ace. '^'^ The 
clean surface acts as a photocatalyst when exposed to 
UV-light. As the intensity of solar radiation on earth 
is maximal for visible light it is important to mod- 
ify the frequency range of the photocatalytic response, 
which is possible by doping with transition metal atoms^ 
and transition metal adatoms/ respectively. The use of 
adatoms is favorable as the 'spatial separation' of the 
excited holes and electrons improves the photocatalytic 
response.^ Adatoms modify the electronic structure, re- 
duce the electronic gap and thus optimize the Ti02 sur- 
face for the use in gas sensors^ and solar cells. ^ In this 
context the Fe-rutile composite is important as it can 
be used for ammonia synthesis. ^^ In addition, the use of 
magnetic transition metal atoms such as Fe integrates 
the spin degree of freedom into the system, which allows 
for new functionalities. For instance, the Fe/Ti02 sys- 
tem has been discussed in the context of diluted magnetic 
semiconductors^^ and a magneto-electrical coupling has 
been detected in Fe-Ti02 systems experimentally.'^ 

General trends of low-temperature growth of (transi- 
tion) metal atoms on the rutile surface are reviewed from 
an experimental point of view in Ref s . [Tl and [T3l and refer- 
ences therein. The growth mode of different metal atoms 
has been classified into two different modes. On the one 
hand, group I and II metal atoms as well as transition 
metal atoms with a small number of valence electrons 
up to Cr interact strongly with the Ti02 surface. Due 
to the large oxidability of these elements, charge from 
these adatoms is transferred to Ti at the surface. As a 
consequence, stable metal-surface bonds form which lead 
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to rather flat metallic overlayers. Also, the formation of 
ternary compounds is likely. On the other hand, if the 
oxidability of the metal atom is weak, e.g., for Cu or Zn, 
no charge is transferred to the surface and the metal- 
substrate interaction is weak. For this reason, no flat 
metal films can be stabilized on the rutile surface. 

The Fe-Ti02 interaction is on the edge between these 
growth modes and the energy of formation of Fe-Ti02 
composite systems is in the same range as the Fe-Fe 
interaction and the formation of Fe clusters. ^^ Indeed, 
different growth modes of Fe on Ti02 have been ob- 
served experimentally, depending on the preparation 
technique used and growth parameters such as temper- 
ature and Fe flux. The formation of fia t Fe ov erlayers 
has been observed at low temperatures IHMEU In con- 
trast, three-dimensional growth and Stranski-Krastanov 
growth (three-dimensional growth after the first layer has 
been formed) have been reported.^^^^ 
Therefore, a more profound understanding of the Fe- 
Ti02 interaction, which can be obtained from ah initio 
simulations, is important in order to optimize the atomic 
and electronic interface structure. Although, the influ- 
ences of temperature and growth kinetics are important, 
simulations at T = K are capable to describe the basic 
mechanisms of the interface formation. So far, system- 
atic theoretical studies on the adsorption of nonmagnetic 
metals such as Mo,^ Pd,^ Rh,!^ and single Cu, Ag and 
Au atoms exist. ^^ However, with respect to iron, only 
the adsorption of single Fe atoms on two different sur- 
face positions ( which are energetically most favorable 
for adsorption of Mo and V) is discussed in Ref. [20| It 
has been predicted in that work that single Fe atoms can 
bind to atoms in different surface positions as long as 
three O neighbors (distance < 3A) and an additional Ti 
bond towards the metallic adatom exist simultaneously. 
The adsorption of Fe-0 molecules has been discussed in 
Ref. 7, It has been highlighted that different adsorption 
sites are thermodynamically stable and that the photo- 
catalytic response is enhanced as the band gap is reduced 
by 0.3 eV because of Fe states near the valence band max- 
imum. 
In the present work, a more detailed study on the adsorp- 



tion of single Fe atoms and a bundle of Fe atoms will be 
presented. The latter setup is essential in order to model 
the growth mode at higher Fe flux or after the deposition 
of the first atoms. We therefore present a first investiga- 
tion of an increasing Fe coverage and of the adsorption 
of small Fe clusters on the Ti02 surface in order to fill 
this gap. It is shown that the Fe-Fe interaction reduces 
the surface wetting for an increasing number of Fe atoms 
deposited at the same time. In addition, the adsorption 
of the first Fe atoms is crucial for the further growth pro- 
cess and the interface symmetry. Thus, our results can 
be used to discuss many facets of experimental results. 
Furthermore, to our knowledge the magnetic properties 
of the rutile surface with adsorped Fe have not been dis- 
cussed on a first-principles level so far. We show that the 
magnetic moments of Fe are stable at the interface. In 
addition, finite magnetic moments are induced at Ti and 
O atoms in the interface via localized hybrid states. Our 
results suggest a similar magneto-electrical coupling as, 
e.g., discussed for Fe-BaTiOs interfaces in Ref. 1211 

The paper is organized as follows. In Section [ll] we 
present the computational methods used in this work 
which is followed by a discussion of the atomic and elec- 
tronic structure of pure rutile and its (llO)-surface in 
Section [nij The adsorption processes of single Fe atoms, 
which corresponds to the experimental setup of low- 
temperature adsorption at low Fe pressure, is discussed 
in Section IV In Section |V| the influence of the Fe-Fe 
interaction on atomic and electronic structure of the Fe- 
Ti02 interface will be discussed for larger Fe coverage 
and for attached clusters. The magnetic properties of 
the material are discussed in Secti on \VJ\ Summary and 
conclusion can be found in Section I VIII 



II. COMPUTATIONAL DETAILS 

The plane wave pseudopotential code VASF^has been 
used for self-consistent electronic structure calculations 
using projector augmented wave potentials^ and the 
generalized gradient approximation of Perdew, Burke, 
and Ernzerhoi"^^. It is a well known fact that the band 
gap of oxidic materials is underestimated by this ap- 
proach. Partly, this failure can be attributed to the un- 
derestimation of the correlation between electrons in lo- 
calized d-states, an error which can be corrected by the 
GGA-\-U approach. For oxygen deficient Ti02, an on- 
site interaction of/7 = 5.5±0.5eV has been determined 
on an ab initio level. ^^ Choosing U = 5 eV (J = eV) 
the influence on the atomic and electronic structure is 
tested by applying the GGA-\-U approach in the formu- 
lation of Dudarev.^^ An energy cutoff of 500 eV and an 
energy convergency of 10~^ eV guarantee highly accurate 
results. In addition, the relaxation of the ions was car- 
ried out until the forces were converged to 0.01 A/eV. 
The /c-mesh has been constructed with the Monkhorst- 
Pack scheme^ and the number of points used has been 
adjusted to the system size, e.g.,17x 17x17 /c-points have 
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FIG. 1. (Color online) (a) Structure of bulk rutile. Relative 
ion size corresponds to covalent radii. Light blue: Ti, dark 
red: O. Arrows mark the two different kind of Ti-0 neighbors: 
shortest equatorial bond Ti-Oeq and apical Ti-Oap bonds. 



been used for bulk rutile whereas the mesh has been re- 
duced to 4x4x1 for the investigation of Fe atom adsorp- 
tion on the rutile surface. A Gaussian-smearing of the 
electronic states of 0.1 and 0.05 eV has been used during 
the optimization of the atomic structure and for the de- 
termination of the electronic structure, respectively. Lo- 
cal magnetic moments were obtained by projecting the 
wave functions onto spherical harmonics within spheres 
of 1.32 A, 0.82 A, and 1.30 A for Ti, O, and Fe atoms, 
respectively. 

In order to reduce the computational effort, a minimal 
number of electrons has been treated as valence including 
the Ti Sd'^s^-, O 25V-, and Fe 4d^45^-states. Although, 
the missing p electrons in the Ti valence base may lead to 
small errors, ^^ the main atomic relaxation as well as the 
most important features of the electronic structure and 
the polar character of Ti02 are well accounted for by 
this approach.^ '^ ' ^^'^^ ' At least 10 A of vacuum was used 
to prevent interactions between the periodically repeated 
slabs in case of free and Fe covered surfaces. 



III. FEATURES OF RUTILE 

The rutile morphology of Ti02 possesses P42/mnm 
symmetry. Beside lattice constant and c/a-ratio of the 
tetragonal distortion, the internal parameter u deter- 
mines the Ti-O distances. Each Ti atom is octahe- 
drally surrounded by O atoms with four nearest equa- 
torial neighbors (Ti-Oeq) and two next-nearest apical 
neighbors (Ti-Oap), see Fig. fl] The material properties 
obtained with the computational tools used here are in 
sufficient agreement with previous experimental and the- 
oretical work, see Tab. |Tj The atomic volume is slightly 
overestimated because of GGA potentials and the well 
known underestimation of the electronic gap is present. 
The description of the gap can be improved by the in- 
troduction of a Hubbard term for the Ti c^-states, but 
simultaneously the description of the structural proper- 
ties worsens. 
The Ti-0 bonds in the pure oxide are mainly ionic but the 



TABLE I. Lattice constants a, c, and internal parameter u, distances of the equatorial (TiOeq) and apical (TiOap) Ti-0 
neighbors in A (cf. Fig. [I]), and static charge per atom Q. We note that the static charge is not a well defined quantity and 
thus different experimental and theoretical approaches yield quantitative different results. 
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formal ionicity of Ti^+02~ is considerably reduced, as Ti 
d- and O p-states hybridize, cf. Ref. [29| The density of 
states (DOS) obtained in Fig. |3] agrees qualitatively with 
previous experimental and theoretical results. ^^ ^^ The 
valence band maximum consists of O p-states, which are 
oriented within the (100) /(OlO) plane, while a significant 
occupation of Ti (i-states with e^-character appears in the 
rest of the valence band. The conduction band is com- 
posed of Ti d-states which possess mainly t2^-character at 
the band minimum. If a Hubbard U-teim is applied, the 
conduction band minimum is shifted to higher energies 
and the hybrid peak which consists of O s-Ti d-states, is 
slightly shifted to lower energies. However, no qualitative 
modification of the atomic and electronic structure has 
been observed if a [/ correction is imposed. 

In the following, the properties of the Ti02 (110)- 
surface, which are essential for the adsorption process, 
will be summarized. We refer to Refs. l29l |3l [371 and|38] 
for further details of the atomic and electronic structure 
of this specific surface. The (llO)-surface is the ener- 
getically most favorable low index surface as it is charge 
neutral and the coordination of the surface atoms is close 
to the bulk values. Within each surface unit, only one Ti 
and one O atom are undercoordinated as rows along the 
[001] -direction of 5-fold and 6-fold coordinated Ti atoms 
alternate (Tisc/Tiec), see Fig. |2] In the remainder of the 
paper the nomenclature Tisc/Tiec will be used for the 
Ti atoms in the surface layer 1 and the corresponding 
atoms in the layers below. Charged TiO+ and bridging 
O^ layers alternate in [110] -direction and each surface 
unit, i.e. each monolayer (ML) consists of two Ob and 
one TiOp layer, see Fig. [2] The surface is slightly buck- 
led, as the undercoordinated Ti^c and bridging O (Ob) 
atoms relax inwards whereas Tiec and O atoms within 
Ti-0 planes (Op) move outwards, see arrows in Fig. l2J 
The Ti-0 distances along the surface normal alternate 
with the surface distance, whereas the modifications of 
the Ti-0 distances within the surface plane are minor. ^^ 
The underlaying mechanism for these modifications of 
the Ti-0 bonds along the surface normal is the reduced 
coordination of Ob and Ti^c atoms in the surface layer. 
For these atoms, the energy difference between Ti d- and 
O p-orbitals, which are aligned along the surface normal, 
is reduced, the hybridization and the bond strength in- 




FIG. 2. (Color online) Atomic structure of a free 3 monolayer 
thick rutile (llO)-film. Arrows illustrate the atomic relaxation 
at the surface. Relative ion size corresponds to covalent radii. 
Ti: light blue, O: dark red. The most favorable adsorption 
position for a single Fe atom is indicated as well. 



crease resulting in reduced Ti-0 distancesJ ^^ * ^^ * ^^ * This 
modification of the bonds is accompanied by a charge 
transfer along the surface normal, see Tab. [Hj For ex- 
ample, the charge of the topmost Ob atom is reduced 
by 14% (17%) within our Bader analyis (the Mulliken 
analysis in Ref. 34) , while the charge of the Ti atoms in 
the first layer increases by 2% (4%). This means that 
local dipole moments develop along the surface normal. 
Since an overall dipole moment along the surface normal 
would destabilize the surface, a stability saving opposite 
charge transfer takes place for the Ti-0 bonds below, 
which results in reduced Ti-0 hybridization and enlarged 
Ti-0 distances for the corresponding Ti-0 bonds. Due to 
the large polarizability of Ti02, the atomic surface relax- 
ations discussed here extend several ML into the surface 
and are thus truncated for thin films, and hence induce 
finite-size effects in the electronic structure. These finite 
size effects open up the possibility of "fine-tuning" of the 
electronic structure of the Ti02 surface by using ultra- 
thin films. In order to separate the general trends of the 
Fe-Ti02 interface from the thickness dependent modula- 
tions, a detailed discussion of the influence of film thick- 
ness on the atomic and electronic structure is necessary. 
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FIG. 3. (Color online) Layer resolved density of states (DOS) 
per Ti02 formula unit, (a) 5 monolayers Ti02 with two fixed 
bottom layers with [7 = 5 eV; (b) as (a) with [7 = eV; 
(c) 3 monolayers with one fixed bottom layer for the case 
[7 = eV. In each case the black dotted line corresponds to 
the bulk reference. Light gray: Bottom layer; Thin blue line: 
central layer; Thick red line: surface layer. 



The atomic relaxations of ultrathin Ti02 films depend on 
whether the films possess a center of inversion (odd num- 
ber of layers) or not (even number of layers). In the first 
case, the amplitude of the discussed oscillation of the in- 
terlayer distances shrinks with increasing film thickness, 
while it increases for an even number of layers. In order 
to reduce the influence of the parity of the number of 
layers, it is convenient to fix the ionic positions at the 
bottom. Here, we restrict ourselves to the discussion of 
an odd number of layers and present a systematic com- 
parison of completely relaxed films and films with fixed 
bottom layers. 

The Ti-Ob distances which we obtained for different film 
thicknesses are listed in Tab.jTTTj For both, free films and 
fixed bottom layers, all Ti-Op distances as well as the 
Ti-Ob distances within the surface are sufficiently con- 
verged with respect to the thickness for 3 ML. Similarly, 
the charge state of the topmost Ob is nearly independent 
of the film thickness, see Table [ill ^^^ ^^^ films, the con- 
vergency of the atomic charges with respect to the film 
thickness is much slower. In this case more than 5 ML 
are required in order to obtain converged charge distribu- 
tions. Although, the overall ionic character of the Ti-0 



TABLE II. Static charge distribution per atom [|e|] within 
the topmost free layers for different film thicknesses for free 
films and films with fixed bottom layers. In the last row the 
tight-binding charge distribution for a free film is listed for 
comparison 



Layer 1 (surface) 
ML Ob Op Tiec Tisc Ob 



t - 



-0.86 -LOl L96 1.96 -L03 

-0.86 -1.02 1.96 1.94 -1.00 

-0.93 -1.06 2.08 2.08 -1.10 

-0.82 -1.06 1.96 1.99 -1.01 

-0.86 -1.02 1.96 1.95 -1.01 



-1.10 -1.17 2.35 2.29 -1.19 



Layer 2 (surface -1) 
Ob Op Tiec Tisc Ob 



-0.98 -1.00 1.98 1.99 -0.95 

-0.98 -1.00 1.98 1.98 -0.96 

-1.05 -1.06 2.11 2.11 -1.04 

-0.95 -0.99 1.92 1.99 -0.95 

-0.99 -1.00 1.97 1.98 -0.96 



2.34 2.37 



^Bulk: 1.99/-1.00;R/ = 5 eV: Bulk: 2.11/-1.06; 

^ :ief[39] Bulk +2.36/-1.18;pFixed bottom layers; pFully relaxed. 



TABLE III. Left: Deviations of the Ti-Oeq distances com- 
pared to a 9 ML film with a fixed bottom layer (%). Right: 
electronic band gap as function of on the number of layers N 
without U correction. 



Distances 






Egap 




Layer 1 Layer 2 










N Tiec Tisc Tiec Tisc 


Total Surface Ref. \TB Ref. EJ 


3 0.5/-1.0 -0.5 0.8 -2.7/3.3 


1.05 


1.05 


- 


- 


c^ 0.5/-1.0 -0.8 1.0 -3.5/4.2 


1.20 


1.20 


1.27 


1.45 


5 0.3/-0.2 -0.3 0.0 -0.5/0.7 


1.1^^ 


1.36 






^^ 0.3/-0.2 -0.3 0.3 -1.2/1.2 


1.43 


1.43 


1.58 


1.71 


7 0.0/ 0.0 0.0 0.0 0.0/ 0.3 


1.19 


1.64 


- 




i' 


[ 0.0/ 0.0 -0.3 0.0 -0.4/ 0.4 


1.51 


1.51 


- 


1.77 


i 


' 0.0/ 0.0 0.0 0.0 -0.1/0.2 


- 


- 


- 


1.85 



ulk: Egap = 1.93 eV;[^ully relaxed; QVith U: Egap =1.54 eV 



bonds is enlarged by a /7 term, the spatial distribution of 
the charge and its convergency with respect to the num- 
ber of layers is practically not modified when U is taken 
into account. 

As the topmost Ti-O bonds are most important for the 
adsorption processes, qualitative results can already be 
expected for ultrathin film of 3 ML thickness. However, 
the Ti-O distances along the surface normal in the center 
of the film show considerable deviations between films of 
3 and 9 ML thickness. As a consequence, the electro- 
static potential of the central Tisc atom is modified for 
the free film which results in a considerable reduction of 
the band gap within the whole film, as Tisc s tates below 
the conduction band edge are induced, D^*^ see Fig. [s] 
Although, the gap opens up with the number of layers, 
it is still underestimated for 7 ML.^ 
If the symmetry of the film is broken by fixing the bot- 
tom atoms, no such overall reduction of the bandgap ap- 
pears. But, artificial surface states are induced for the 
fixed atoms. Most notably, a Tisc-state below the con- 
duction band minimum appears in the fixed layer which 
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FIG. 4. (Color online) Change in total energy for adsorption 
of Fe on the rutile (llO)-surface at positions 1-7 relative to 
the ground state (position 1), cf. Fig. [2] 



has mainly (ia,2_ ^2 -character and hybridizes with the Tisc 
atom above. Furthermore, a Ti-0 hybrid peak appears 
at the valence band maximum. The energetically highest 
occupied states within this peak are p-states of the fixed 
Ob atom. The overall gap is therefore smaller than for 
free films. However, the surface states decay rapidly with 
the distance towards the fixed surface and the conver- 
gency of the electronic structure in the free surface with 
respect to the film thickness is considerably improved. 
Only for ultrathin films of 3 ML thickness, hybridization 
between these artificial surface states and the states of 
the central layer appears, see Fig. Is] (c). Thus, small 
deviations in the Fe adsorption process between thick 
films and a film of 3 ML thickness are possible. Qualita- 
tively, the discussed modification of the electronic struc- 
ture with the number of layers is not modified if a Hub- 
bard term is imposed on the Ti d-states. 

In summary, we could show that the atomic arrange- 
ment and the charge distribution within the uppermost 
free layer have practically converged for a film of 3 ML 
thickness which is therefore sufficient for the study of ba- 
sic adsorption mechanisms. In addition, the convergence 
of atomic and electronic structures can be improved as 
the bottom layers are fixed. 



IV. ADSORPTION OF A SUB-MONOLAYER OF 

IRON 

In this section we describe the modeling of the adsorp- 
tion process of Fe atoms at low-temperatures and low Fe 
flux. Under such experimental conditions individual Fe 
atoms are adsorbed on the surface and Fe-Fe interactions 
as well as diffusion of Fe atoms on or into the surface can 
be neglected. This process is modeled at zero tempera- 
ture by posing single Fe atoms on the seven different ad- 
sorption positions of the rutile surface (cf. Figs. [2] and [4| 
and by relaxing the atomic positions under the constraint 
of a fixed bottom layer and fixed in-plane positions of Fe. 
In order to minimize the interaction between periodic Fe 



TABLE IV. Characterization of adsorption positions 1-7 in 
Fig. [4] Listed are energy differences A^ (eV/supercell) rel- 
ative to position 1, the change of atomic charges in the Fe 
neighborhood AQ (|e|), nearest neighbor Fe distances d (A), 
and magnetic moments /x (/iB/atom). 



6 



7' 



AE 



0.00 0.37 0.75 1.18 
0.00" - 0.72"^ - 



1.90 1.94 3.21 



AQ 



Fe 


1.22 


1.10 


1.19 


1.02 


1.16" 


- 


- 


- 




0.85^ 


- 


- 


- 


Tiec 


-0.08 


- 


- 


- 


Ti5c 


- 


- 


-0.05 


-0.07' 


Ob 


-0.14^ 


-0.15" 


-0.12 


-0.14 


Op 


-0.14 


- 


-0.09" 


-0.05 



0.81 0.89 0.55 



o.7r 



0.13 



-0.12" 



-0.05" 



0.46' 



-0.10 



Tiec 3.20 2.72 2.91" 2.67 3.50" 

Tisc - - 2.75 2.8" - 2.57" 2.63 



Ob 2.50" 1.80" 1.78 2.07" 



1.73 
3.4" 



3.66" 3.5" 



Op 3.1 3.2" 1.99" 



1.93 



3.9" 1.92" 3.0" 



Fe 


3.3 


3.3 


3.3 


3.1 


3.1 


3.0 


3.6 


Tiec 


-0.1 
0.5" 


-0.01 
0.4" 


0.01 


-0.02 


0.01 


- 


0.01" 


Tisc 


- 


- 


-0.04 


-0.18 


_/ 


-0.13 


-0.06 


Op 


0.08 


- 


0.06 


0.04 


- 


0.06 


0.02" 


Ob 


0.13" 


0.15" 


- 


0.06 


0.08 


- 


- 



"Ref. 20; ^4 ML; "two bonds of equal length; 

^0.05 electrons are transferred to Tiec in the [001]-row; 

(d(Ti-Fe)=4.9 A); ^four bonds of equal length; 

•^0.02/iB are induced at all Tisc atoms within the first layer. 



images, a supercell of (2 x 3) surface units is used which 
ensures a distance of 13.2 A (8.9 A) along [001] ([110]) 
between the periodic Fe images. The computational ef- 
fort is minimized by using a minimal film thickness of 3 
ML with fixed atomic positions in the bottom layer. In 
a similar manner the adsorption of single atoms on two 
different surface positions (which are energetically most 
favorable for Mo and V) has been discussed in Ref. [20l 
The detailed energy landscape which has been obtained 
in the present investigation is shown in Fig. [4j An Fe 
atom attached between two Ob surface atoms and one 
Op surface atom (position 1) is most favorable. In this 
configuration, Fe has tetrahedral O environment which 
also occurs in various oxides. 

The adsorption is 0.37 eV less favorable in case that Fe 
is placed above the Tiec atom, with two Op neighbors 
(position 2). For position 3 the Fe atom has two Op 
neighbors. This configuration is 0.75 eV higher in energy 
than position 1. All other configurations are more than 
1 eV higher in energy; whereby the most unfavorable po- 
sition is on top of the Tisc atom, see Tab. IV and Fig. [4] 




■ Clean surface 

■ Site 1 
H Site 2 
- Site 3 



FIG. 5. (Color online) Modification of the density of states 
for different adsorption positions, (a) Fe at different adsorp- 
tion sites; (b) corresponding Ti02 surfaces. Thin vertical 
lines mark Fe-Ti-0 hybrid states for each configuration. The 
DOS of the free surface is shown as reference, which has been 
aligned with the lower valence band edge. 



Similar energy differences between position 1 and 3 have 
been obtained in Ref. 20 for a rutile film of 9 ML. This 
confirms the fast convergence of adsorption energy with 
film thickness which has been predicted for the case of 
Pd adsorption, see Ref. |18| 

For all adsorption positions, the amount of charge (AQ) 
is transferred from Fe to the substrate, in agreement with 
photoemission measurements of Diebold et al. where 
Fe^+ has been detected on top of Ti02.^^ This charge 
transfer is correlated with the adsorption energy, see 



Tab. IV The largest AQ arises for position 1, while AQ 
is a factor of two smaller for the least favorable site 7. 
This correlation indicates the formation of ionic bonds 
between Fe and the surface, in agreement with the gen- 
eral trends of transition metal adsorption on Ti02.^^ 
However, no linear dependence between charge transfer 
and adsorption energy exists, e.g., AQ of position 3 is 
similar to that of position 1 despite the energy differ- 
ence. Most notably, different atomic relaxation energies 
due to the induced distortion of the Ti02 structure dis- 
turb the linear dependence. For example, only Ob atoms 
next to Fe relax out of the surface for position 1 and thus 
the energy costs for the surface relaxation are minor. In 
contrast, a larger energy penalty for the surface relax- 
ation has to be expected for site 3 as Ob and Op atoms 
relax towards Fe. Similarly, the atomic rearrangements 
of position 5 are restricted to the Ob atoms next to Fe, 
whereas the Tisc and Op atoms relax relative to each 
other for position 6. 

The adsorption mechanism can be understood by a closer 
look at the involved charge transfer and electronic states. 
The charge state of Ob calculated by Bader partitioning 
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FIG. 6. (color online) (a)-(b) The influence of Hubbard corre- 
lations on the density of states for adsorption sites 1 (Fig. E| 
and a film thickness of 5 ML. (a) Fe and (b) Ti02. Metal 
induced gap states (MIG) I and II are marked. The vertical 
dotted line marks the employed separation between free and 
bound charge, see text. 

(c) Schematic sketch of the photoemisson spectrum adopted 
from Ref. |40l The spectrum of the clean surface (black) is op- 
posed to a reduced surface (purple, dashed) and an increasing 
Fe coverage (shades of blue) . 



is reduced by 0.14 |e| per atom at the clean rutile sur- 
face. As a consequence, the electronegativity of these 
Ob atoms is enhanced. This promotes the formation of 
ionic Fe-Ob bonds accompanied by charge transfer to the 
film. Adsorption of Fe atoms next to Ob is therefore most 
favorable whereas position 7 is least favorable, since all 
Fe-0 distances are larger than 3 A which prevents the for- 
mation of strong bonds. In the film, the injected charge 
spreads towards neighboring Ti atoms, see Tab. |lVl In 
agreement, a reduction of the formal Ti^+ state to Ti^+ 
has been observed experimentally.^^ 
Figure [5] illustrates the modification of the DOS for ad- 
sorption positions 1-3. Within the entire valence band 
Ti-Fe and 0-Fe hybrid states form. The states in the 
lower part of the valence band are marked by solid verti- 
cal lines. For positions 1 and 3 Fe-0 hybrid states form 
at the bottom of the valence band and Ti-Feec, Ti-Fesc 
hybrid states at ~-6.2 eV appear, respectively. Similarly, 
the next-nearest Feec atoms hybridize with Fe states re- 
sulting in a hybrid peak at —6.5 eV for site 2. Besides, 
metal induced gap states (MIG) form in the gap of the 
rutile surface (states I below the conduction band and 
states II above the valence band as marked by dotted 
vertical lines in Fig. Isl in full agreement with experi- 
mental results, see Fign6](c). 

For sites 1 and 3 the Fe state I hybridizes with neigh- 
boring Tigc and Tisc atoms, respectively and a state is 
induced below the conduction band edge. This state 
can be assigned to a Ti^+ surface state as obtained in 
experimentaP^ and theoreticaP^ investigations of the re- 
duced rutile surface, see Fig. |6](c). 

For Fe on adsorption site 2, the charge transfer towards 
Ti is much smaller and no Ti defect state forms which 
could hybridize with the Fe state. The differences in the 



Ti-Fe interaction and the corresponding charge transfer 
can be traced back to the atomic arrangement. For sites 
1 and 3 the neighboring Ob and Op atoms mediate an 
indirect couphng between Ti and Fe and promote a large 
charge transfer. In contrast, Fe is positioned directly 
above Tiec at site 2, which enforces a direct Fe-Ti inter- 
action without O mediated interaction. Such Fe-Ti inter- 
action seems to be unfavorable and the charge transfer to 
the surface Ti atom is hindered. Analogously, no direct 
Fe-Tisc interaction appears if the Fe atom is positioned 
above Tisc (site 7). 

A second MIG (state II) is visible at the valence band 
edge for all three adsorption sites, which is mainly of 
O p-type and localized at the Ob neighbors, see Fig. |5| 
Thus the state observed experimentally in the same en- 
ergy window in photoemission spectra^^ is most likely 
an Fe-0 hybrid peak. With increasing adsorption en- 
ergy, the energy of this Fe-Ti^+ state increases and ap- 
proaches the valence band edge for the less favorable po- 
sitions. Thus, the reduction of the band gap depends on 
the strength of the interface bonds. 

In summary, the interpretation of experimental results in 
Ref. 40 that Fe bond charges transfer via an 0-mediated 
charge transfer to interface Ti atoms, could be confirmed. 
For all adsorption sites, approximately 60% of the in- 
jected Fe charge is transferred to its neighboring atoms. 
In addition, free carriers occupy the lowest Ti states 
throughout the Ti02 film and the valence states are 
shifted to lower energies, see Fig. [6] This charge spillout 
is a common failure of standard density functional the- 
ory potentials as the underestimated band gap enforces 
a wrong alignment of the Fermi energy at the metal- 
insulator interface. "^^ Besides, no major deviations of the 
upper valence states are enforced in comparison to the 
free surface, see Fig. [Sj Particularly, the gap below these 
spurious free carriers is conserved and the free and bound 
charges are well separated, e.g., at ~ —0.4 eV for site 1, 
see Fig. ^'^^ 

As the g ap depends on the thickness of the Ti02 film 
(c/. Tab. |III[ ) also the magnitude of the charge spill-out 
differs, compare Fig. [5] and Fig. [6] However, the wrong 
alignment cannot be prevented even if an even number 
of Ti02 layers is used for which the gap is slightly larger 
than for bulk. 

A more accurate modeling of the interface would require 
an improved approximation of the exchange correlation 
functional such as the hybrid functional Bl-WC,^^ which 
is out of the scope of the present work. The use of a 
Hubbard term allows for a qualitative cross-check of the 
obtained results without the spurious band alignment. 
Indeed, the charge spill-out is reduced if [/ = 5 eV is 
applied to the Ti d-states previous to the optimization of 
the atomic structure while the main features of the DOS 
are not modified, see Fig. [6] A more detailed discussion 
on the modification of the DOS is not convenient as the 
value of U has been determined for the deviating setup of 
clean reduced surfaces and is not necessarily transferable 
to the present set up .^ 
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FIG. 7. (Color online) (a)- (b) Layer resolved density of states 
for adsorption site 1 and different Fe coverages. A Hubbard 
term of [/ = 5 eV has been applied to the states marked by 
dotted lines, (a) Fe; (b) Ti02; blue shaded region: MIGs of 
type I and II (c) topview onto the interface. 
* States shifted by 0.2 eV to lower energies for better comparison. 



The hybrid states which have been discussed so far are 
localized at the interface and are not related to the spuri- 
ous charge spill-out. Most notably, the description of the 
MIGs is in good agreement with experimental results, see 
Fig.|6] First of all, interface states I appear at ~ —0.8 eV 
and '^ — 1 eV below the Fermi level within the calculated 
DOS and the photoemission spectrum, respectively. Sec- 
ond, the interface states II are located directly above the 
valence band edge in both cases. 

In summary, the wrong band alignment and the related 
charge spill-out may modify the details of the atomic and 
electronic structure at the interface. However, the main 
qualitative trends are not affected. 



INCREASING COVERAGE OF RUTILE (110) 
WITH IRON ATOMS 



So far, the adsorption of single Fe atoms has been 
modeled which corresponds to a coverage of 0.1 ML of 
bulk Fe. The Fe-Fe distances employed suppress Fe-Fe 
interactions. Besides charge spill-out and localized MIGs 
(Fig.[6|, the insulating nature of the Ti02 surface is con- 
served. If the coverage increases to 0.4 ML, i.e., one Fe 
per (1x1) surface unit, the Fe-Fe distances shrink to 
2.9 A and 6.6 A along [001] and [IlO], respectively. The 
atomic and electronic structure for this coverage is illus- 
trated for site 1 in Fig.JT] Since the nearest Fe-Fe distance 
is approaching its bulk value in this case, Fe states hy- 
bridize with each other. The DOS peaks broaden and 
interface states I and II are spread over the main part 
of the gap, see the shaded region in Fig. [7[ As for the 
coverage of 0.1 ML, the states can be classified as Ob and 
Op-Fe hybrid states above the valence band edge (MIGs 
I) and Tiec-Fe hybrid states below the conduction band 
edge (MIGs II). All states in the shaded regions are lo- 




TABLE V. Distances (A) of Fe towards its nearest neighbors 
(d < 3 A) for the interface configurations in Fig.pland relative 
energies in eV/(Fe atom). 
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FIG. 8. (Color online) (a)- (c) Different adsorption scenarios 
for 0.8 ML Fe on rutile (110) (d) Density of states of Fe on 
rutile (110) for 0.4 ML Fe (position 1) and 0.8 ML Fe (con- 
figuration (a)) in comparison to bulk Fe. 



calized in space and are qualitatively not modified by a 
Hubbard term. Analogously, increasing weight and width 
of the interface states with increasing Fe coverage have 
been detected experimentally/^ ^^ see Fig. p^(c). 
The weight of the lower third of the valence band in- 
creases while states of the upper third are depopulated for 
increasing Fe coverage in agreement with experiment.'^ 
For example, Fe-Ob states in the shaded region around 
—7 eV gain weight and a Fe-Tiec hybrid state builds up 
at ~ —4.4 eV, see Fig.[7|(a)-(b). The system gains energy 
by the formation or strengthening of hybrid states due to 
the adsorption of additional Fe and, hence, all states are 
shifted by ~ 0.3 eV to lower energies. 
By this mechanism, an increasing fraction of the lowest 
Tisc states at the conduction band edge are occupied. 
Since the states in the central Ti02 layer are not affected 
by this modified band alignment, see Fig.[7|(b), no quali- 
tative modifications of the interface properties due to the 
enhanced spill-out have to be expected. 
The further increase of coverage to 0.8 ML, i.e., two Fe 
atoms per surface unit, corresponds to a closed layer 
which is strained to the lattice constant of Ti02 . In order 
to investigate this coverage, the crystallographic orienta- 
tion TiO2(110)[110]/Fe(001)[010] is assumed which has 
been determined for Fe nanocrystals in a rutile matrix.^^ 
Three representative local minima of the energy land- 
scape are shown in Fig. [s]^^ A homogeneous distribution 
of Fe is not favorable because of tensile strain of 16% 
acting along [110] and because Fe atoms reduce their 
[110] -distances. This relaxation is most pronounced for 
configuration (c) for which both adatoms have initially 
been deposited next to site 3 (Fig. H). 
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After the formation of Fe-Ob and Fe-Fe bonds, the atoms 
are trapped in local minima of the energy landscape 
as all further Fe-0 and Fe-Ti distances are larger than 
3 A. The Fe-Fe distance along [110] is the same range as 
for the free dimer for which distances between 2.01 and 
2.26 A depending on the spin configuration are given in 
literature. ^^ 

For configuration (b) the Fe atoms have initially been 
put at positions 2 and 6. An increasing amount of Fe-0 
and Fe-Ti bonds can form which lowers the total energy, 
although energy penalties have to be expected for the in- 
creasing Fe-Fe distance and the enhanced relaxation of 
the Ti02 surface in comparison to (a). 
If the Fe atoms are initially placed far from Ob atoms the 
adsorbed atoms are not trapped at the surface during the 
first relaxation steps and can thus arrange in a more fa- 
vorable relative alignment, i.e., one additional Fe bond 
per surface unit is formed and each Fe atom possesses 
four Fe neighbors in comparison to 3 Fe-Fe bonds for con- 
figuration (b) and (c) , which results in a large energy gain 
of 0.5 eV. The Fe positions correspond approximately to 
positions 3 and 6 after relaxation. In comparison to the 
adsorption of single Fe atoms, the Fe-0 bonds are weaker 
due to the new Fe-Fe interaction i.e., the Fe-0 distances 
increase, compare Tab. IV and[Vl 

Based on experimental results the formation of metal- 
lic Fe has been predicted for a coverage of 0.7 ML^^^ 
and 3.0 ML,^^ respectively. Although, Fe states start to 
form bulk like bands for a coverage of 0.8 ML, the pro- 
nounced Fermi edge of metallic Fe has not yet formed, 
see Fig. [I (d). 

Further increase of coverage is modeled for a layer-wise 
growth as was predicted in J^, For this purpose, 3 and 5 
strained layers are deposited above the surface which cor- 
responds to 2.4 and 4 ML of bulk Fe, respectively. The in- 
terface configuration shown in Fig. [81(c) has been used as 
starting point for both coverages and the structures ob- 
tained by relaxation are shown in Figure [9] (a) -(b). Most 
notably, only two Ob-Fe bonds per surface unit build up 
which possess a similar bond length as for 0.8 ML cover- 
age. This means that an energy barrier exists for further 
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FIG. 9. (Color online) Atomic relaxation pattern for (a) 
2.4 ML; (b) 4 ML; (c) 1 ML; (d) 3 ML Fe on rutile (110). 
Representative Fe-0 and Fe-Fe distances are included, (a)- 
(b) Fe is strained to the lattice constant of rutile; (c)-(d) Fe 
superstructure, see text. 



FIG. 10. (Color online) Topview on the initial configurations 
of a cuboid Fei4 cluster attached to a (2 x 3) surface cell of 
rutile (110). Atoms within the direct surface are shown. The 
energy differences of the optimized structures relative to (a) 
are given in eV/(Fe atom). 



In summary, for an increasing coverage the interface 
structure depends on the adsorption of the first deposited 
atoms which shows the influence of different kinds of 
growth parameters on the interfaces in experiments. The 
formation of flat films with bcc- morphology is possible ei- 
ther by formation of strained films or by different lattice 
periods along [110] in rutile and Fe. 



V.l. Adsorption of Fe clusters on rutile (110) 



interface bonds, if the first Fe layer is adsorbed in the 
metastable state which is shown in Fig. Is] (c) and which 
cannot be overcome optimizing the structure at T = K. 
The tensile strain for three adsorbed layers in combina- 
tion with the reduced coordination of the surface atoms 
leads to an inwards relaxation of the uppermost Fe layer 
and two buckled Fe layers form instead, compare Fig- 
ure |9Fa). By this relaxation, each Fe atom obtains ei- 
ther three or four nearest Fe neighbors with a distance 
of d ^ 2.3-2.6 A, i.e., in the range of the bulk distance 
(2.45 A). In addition, up to four next-nearest Fe-Fe bonds 
with d - 2.7-2.8 A (bulk value 2.83 A) form. If the film 
thickness increases, the relaxation of the Fe layers is less 
pronounced as the ratio between surface atoms and bulk- 
like coordinated atoms in the center of the film decreases 
and a rather flat film is stabilized, see Figure lojb). 
Besides the formation of strained films, one may think 
of films which are not congruent with the underlying lat- 
tice constant along [110]. For example, a superstructure 
with a ratio of 2:5 between the Ti02 and Fe lattice con- 
stant would reduce the epitaxial strain and the mean Fe 
volume is compressed by only 1.3%.^^ The correspond- 
ing relaxed interfaces for 1 and 3 ML coverage are shown 
in Fig. [9Fc-d). In both cases, the adatoms have initially 
been placed above Tiec and next to Op atoms in order 
to circumvent the unfavorable trapping near Ob atoms. 
During relaxation, the Fe-Fe distances shrink and at least 
one Fe-0 bond is formed for each Fe with bond lengths 
in the range of 2 A for a coverage of 1 ML. The Fe-0 
interaction is reduced, i.e. the Fe-0 distances increase 
if the Fe-Fe coordination is increased by a coverage of 
3 ML. Nevertheless, flat films with good surface wetting 
are obtained. 



Complementary to the adsorption of single atoms, the 
case of clusters attached to the rutile surface is modeled 
in the following. For this purpose, the adsorption of Feis 
and Fei4 clusters on the clean Ti02 surface has been in- 
vestigated. In both cases, (2 x 3) surface units are used, 
in order to reduce the interaction between neighboring 
clusters as schematically shown in Fig. [lO] 
The ground state of free Feis clusters is a Jahn- Teller 
distorted icosahedron with FM alignment of all magnetic 
moments. ^^ If such a relaxed cluster is placed above the 
Ti02 surface with its vertex pointing towards the rutile 
surface, two adsorption trends can be distinguished. If 
the vertex is initially positioned above the Op-Tisc region 
of the surface, the cluster aligns itself in between two Ob 
rows and thus two strong Fe-Ob bonds form at each side 



11 ^a). For the 



of the cluster {d ^ 1.9 A) as shown Fig. 
second setup, the vertex of the cluster is initially placed 
near Ob atoms, see Fig. ITTlb). In this case, two Fe-Ob 
bonds with a length of d ^ 1.9 A build up. In this case 
the system is trapped in a local energy minimum as has 
been discussed for the adsorption of films and a smaller 
number of interface bonds forms, e.g., two (one) Fcvertex- 
Tisc {d ^ 2.4 A) and Fcvertex-Op {d ^ 2.0 A) bonds form 
for (a) and (b) in Fig. [TT] respectively. Therefore, con- 
figuration (b) is 0.3 eV/(Fe atom) less favorable. 
For both configurations, the relaxation of the Ti02 sur- 
face is minor and restricted to direct neighbors of Fe. 
Only the relaxation of one Op and Ob atom towards 
the Fe cluster for (a) and (b), respectively, is significant. 
For cluster (b) the whole cluster tilts towards the sur- 
face allowing for one Fcgide-Op {d ~ 2.0 A), and two 
Fcside-Ob bonds {d ~ 2.0 and 1.96 A).^^ Analogously, 
one Fcside-Tisc bond {d ^ 1.9 A), and one Fcgide-Op bond 
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FIG. 12. (Color online) Density of states of (a) Fe and (b) 
Ti02 of the adsorbed Fei4 cluster (Fig. flTld)). For compari- 
son results for a flat film (Fig.^a)) and bulk Fe are given. 



FIG. 11. (Color online) Adsorption of Fe clusters on rutile 
(110). (a)-(b) Icosahedral Feis (c)-(g) Cuboid Fei4; (c) Initial 
configuration corresponding to (d); (d)-(g) relaxed configura- 
tions sorted by relative energies (eV/ (Fe atom)); magnetic 
moments (/is /atom) are given for Fe if the value differs from 
2.8-3.0 /iB and for representative Ti and O atoms. 



{d ^ 1.9 A) form in case of configuration (a). 
In both cases, the icosahedral morphology of the cluster 
is conserved without a tendency for surface wetting. As 
the closed-shell Feia icosahedron is particularly stable, a 
large energy barrier for the fracturing of such cluster has 
to be expected, which cannot be overcome in the simu- 
lations at zero temperature. 

In order to test whether this energy barrier or the intrin- 
sic strength of Fe-Fe and Fe-Ti02 bonds hampers surface 
wetting, the adsorption of cubic Fei4 clusters with bcc 
morphology has been investigated. Free clusters of this 
geometry are not favorable because of the reduced coor- 
dination of Fe atoms in comparison with an icosahedral 
geometry and therefore the energy barrier for the redis- 
tribution of Fe may be overcome in T = K simulations. 
As in the previous section, the crystallographic orienta- 
tion TiO2(110)[110]/Fe(001) [010] has been assumed. The 
cluster is strained to adopt the lattice constant of rutile, 
and several adsorption positions have been sampled as 
shown in Fig. 
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The deposisted cluster is indeed not 
stable against rearrangement of Fe atoms, see Fig. 11 ^d- 
g). For all tried adsorption positions, the atomic relax- 
ation shows similar behavior. First, Fe-0 bonds form 
with a subsequent collapse of the clusters which results in Fermi edge exists in agreement with experiment. ^=^ 



different distorted clusters with local icosahedron geom- 
etry. For configurations (f) and (g), the cluster is again 
trapped in the local energy minimum next to the Ob 
atoms. Furthermore, the simulation cells for these con- 
figurations possess mirror symmetry, which is conserved 
during the relaxation. This leads to a minor Fe-substrate 
interaction. 

Although, the used conjugate gradient method does not 
reproduce the time-evolution of the system, the unique 
relaxation trends underline that Fe-0 interaction is the 
driving force for Fe adsorption. 

In particular, the low wetting tendency for an unfavor- 
able cluster geometry shows that the formation of dense 
clusters is energetically more favorable than the forma- 
tion of further Fe-Ti or Fe-0 bonds. Thus, a three- 
dimensional growth mode is likely, if clusters instead of 
single Fe atoms are deposited onto the surface. 
The DOS of the adsorbed Fei4 cluster agrees qualita- 
tively with the electronic structure discussed in Sec. |IV[ 
see Fig. 12 On the one hand, the lowest conduction 
band states are occupied in the entire film because of the 
charge spill-out. On the other hand, localized MIGs are 
induced in the shaded energy interval, see Fig. 12 'b). In 
comparison with the adsorption of single atoms, these 
interface states are broadened to a "defect" band which 
fills the gap and the energy intervals of MIGs I and II 
overlap. 

Due to the compact cluster which is formed by relaxation, 
the Fe DOS is more similar to the bulk DOS as in case 
of flat films and bulk like bands form, especially for the 
central atom of the Fe cluster. However, no pronounced 
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FIG. 13. Spin-polarization of Fe in the first Fe adlayer 
(Fig. [9|b)) and in a deposited Fei4 cluster (Fig. 11 'b)) on 
rutile (110) in comparison to the spin-polarization of bulk Fe. 



VI. MAGNETIC PROPERTIES AND 

MAGNETO-ELECTRICAL COUPLING OF 

RUTILE (110)/ IRON 

The adsorption of Fe on the rutile surface incorporates 
the spin degree of freedom and new functionalities may 



be thought of. Table IV lists the magnetic moments for 
single adsorbed atoms. For all adsorption sites the mag- 
netic moment is reduced compared to the value of the free 
atom (4 /ie)- For positions 1-6, the Fe-Ti02 interaction 
reduces the Fe moments by ~ 20% whereas the weak 
bonding in case of position 7 results in a reduction of 
only 10%. Most notably, all Fe moments are larger than 
the bulk value of of 2.2 //b- If the Fe coverage increases 
to 0.4 ML, 0.8 ML, 2.4 ML, and 4 ML in a layer-wise 
growth mode, mean magnetic moments of 3.0, 2.9, 2.8, 
and 2.7 /iB are obtained in the simulations, i.e. the mag- 
netic moments slowly converge towards the bulk value. ^^ 
For adsorbed Feia and Fei4 clusters, the interface atoms 
without strong interface bonds possess magnetic mo- 
ments of 2.9 /iB similar to the moments of 3.1 and 2.8 
obtained for the free Feia cluster in its high- and low spin 
state, respectively.^^ The moments of the central atoms 
are reduced with respect to the high-spin phase of the free 
cluster (2.8 /iB^^) and the low-spin state with an antifer- 
romagnetic alignment of the central atom is not stable for 
adsorbed Fei4 clusters. If strong Fe-0 and Fe-Tisc bonds 
form the corresponding magnetic moments are reduced 
whereas the moments are not modified in case of weak 
interface bonds, see Fig. 11 An exception is the ideal 



icosahedron (g) for which the large oxygen coordination 
of the vertex atom results in an ant i- ferromagnetic align- 
ment of this atom. 

In summary, the magnetic moment of Fe is stable in the 
vicinity of the Ti02 surface which is essential for spin- 
tronics applications. 

Experimentally, a magneto-electrical (ME) coupling has 
been determined for Ti02-Fe composites. ^^ Such cou- 
pling can be mediated by interface hybridization^ 
or spin-accumulation.^ The discussed pathologic band 
alignment and the resulting charge spill-out forbid a 
quantitative estimation of the coupling in the present 



investigation.^^ However, a spin-polarization at the Fermi 
level of 57% is obtained in bulk Fe,^^ which is most likely 
sufficient for a spin-accumulation during the screening of 
the interface charge which would appear under an ap- 



plied electrical field. Figure [13] illustrates the calculated 
spin-polarization of the composite material. At the Fermi 
level, the polarization is reversed with respect to the bulk 
material and a large polarization of -90% is obtained, 
which would allow for spin accumulation at the inter- 
face. 

Complementary to the expected charge transfer and spin- 
accumulation, a significant ME coupling can be expected 
because of the magnetic properties of the Fe-Ti02 hy- 
brid states analogously to the discussion of BaTiOa/Fe 
in Ref. 21, Besides the modification of the Fe moments 
at the interface, magnetic moments are induced in the 
Ti02 film as both the spurious free charges and the MIGs 
at the interface are spin-polarized, compare Figs. [6][7| and 
Tab. IV As the lowest Ti^c <i-states align with the minor- 
ity Fe states, small magnetic moments are thus induced 
at the Tisc atoms in the whole film, which are antiparal- 
lel to the Fe moments. 

For 0.1 ML Fe on adsorption position 1 the magnetic 
properties of the bound charges can be separated from 
these spurious charges by integrating the spin-polarized 
DOS up to —0.4 eV, compare Fig. [5] Within this energy 
interval, only magnetic moments in the direct interface 
exceed the accuracy of 0.01 /iB and the induced moments 
can be traced back to the hybridization with MIGs I and 
II. Magnetic moments of ~ 0.1 fiB are induced at Ob and 
Op atoms for which a Fe-0 hybridization with the Ti 
peak I takes place. The Ti^c atom with a Fe-Ti distance 
of 3.2 A and the next-nearest Tiec atoms hybridize with 
the minority state II and magnetic moments of —0.1 /iB 
and —0.05 /iB are induced, respectively. 
Both, the Fe moments at the interface and the induced 
Ti moments are sensitive to the interface configuration 
and are thus most likely modified if polar displacements 
are induced in the film. For instance, no hybridization 
between the Fe state II and Ti appears for position 2 and 
thus no magnetic moments are induced at the interface 
Ti atoms, see Fig. [5] 

For a coverage of 0.4 ML, states are induced within the 
whole band gap, see Fig. [Tj However, the separation 
of the magnetization up to —0.4 eV is still a good ap- 
proximation for position 1, as the surface states in this 
energy interval decay rapidly with the interface distance. 
In this energy interval magnetic moments of 0.16,— 0.25, 
and 0.07 /iB are induced at Ob, Ti^c, and Op atoms, 
respectively. The induced moments in the Ti02 film are 
increased as more polarized interface states are occupied. 



cf. the induced moments for 0.1 ML coverage in Tab. IV 



In summary, magnetic moments are induced by the inter- 
face bonds, which are sensitive to the Fe-Ti and Fe-0 dis- 
tances and thus a similar ME-coupling as in BaTiOs/Fe 
can be expected. For a more quantitative analysis of the 
ME-coupling strength, accurate hybrid functionals and 
finite electrical fields must be considered in further cal- 
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culations. 



VII. CONCLUSIONS AND OUTLOOK 

We have investigated the adsorption of Fe on the ru- 
tile (llO)-surface by means of ab inito calculations. First, 
the properties of the free rutile (llO)-surface have been 
reviewed. We could show that the slow convergence of 
the surface properties with the film thickness can be im- 
proved if the bottom layers are fixed to their bulk po- 
sition. Already for a film thickness of three layers the 
atomic and electronic structure of the direct surface is 
mainly converged. 

Our investigation of the adsorption of single Fe atoms 
shows that mainly ionic Fe-interface bonds form. The ad- 
sorption next to the topmost bridging O atoms is found 
most favorable as this atom has an enhanced electroneg- 
ativity at the free surface. The Fe-0 bonds mediate a 
charge transfer to the neighboring surface Ti atoms which 
are partly reduced. By the formation of Fe-Ti and Fe-0 
hybrid states, the band gap is reduced which is favorable 
for phot ocathaly tic applications. 

The comparative study of the interface structure for dif- 
ferent start configurations and different amount of de- 
posited Fe confirms the wide spread of experimentally 
found interface structures. On the one hand, a large 
interface coupling and good surface wetting can be ob- 
tained if single Fe atoms are deposited. For this setup 
we found short Fe-interface bonds with a large charge 
transfer from the Fe atom into the surface. On the other 



hand, Fe-Fe bonds reduce the interface coupling for de- 
posited clusters resulting in a low surface wetting and a 
three-dimensional interface structure. Furthermore the 
interface structure depends on the adsorption of the first 
atoms which reflects the influence of the experimental 
growth process. 

Besides the structural and electronic properties, we also 
discussed the magnetic properties of the interface. We 
found stable magnetic Fe moments and a large spin po- 
larization of the adsorped Fe atoms. In addition magnetic 
moments are induced at the Ti and O atoms in the in- 
terface. Most notably, the magnetic interface moments 
depend on the actual interface geometry and thus a hy- 
bridization mediated magneto-electrical coupling is very 
likely. 

However, more sophisticated approximations for the ex- 
change and correlation potential such as hybrid function- 
als have to be used in the future in order to obtain the 
correct electronic gap. In this case a quant it ive estima- 
tion of the photocathalytic properties and the magneto- 
electric coupling will be possible. 
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